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Our attention focuses on the stochastic dynamo equation with non-normal operator that gives 
an insight into the role of stochastics and non-normality in galactic magnetic field generation. 
The main point of this Letter is a discussion of the generation of a large-scale magnetic field that 
cannot be explained by traditional linear eigenvalue analysis. We present a simple stochastic 
| model for the thin-disk axisymmetric aVL dynamo involving three factors: (a) non-normality 

generated by differential rotation, (b) nonlinearity reflecting how the magnetic field affects the 
turbulent dynamo coefficients, and (c) stochastic perturbations. We show that even for the 
subcritical case, there are three possible mechanisms for the generation of magnetic field. The 
first mechanism is a deterministic one that describes an interplay between transient growth 
and nonlinear saturation of the turbulent a— effect and diffusivity. It turns out that the trivial 
state is nonlinearly unstable to small but finite initial perturbations. The second and third 
are stochastic mechanisms that account for the interaction of non-normal effect generated by 
differential rotation with random additive and multiplicative fluctuations. In particular, we 
show that in the subcritical case the average magnetic energy can grow exponentially with time 
due to the multiplicative noise associated with the a— effect. 
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The generation and maintenance of large scale magnetic fields in stars and galaxies has at- 
tracted enormous attention in past years EP-[I] (see also a recent review (Hj)- The main candidate 
to explain the process of conversion of the kinetic energy of turbulent flow into magnetic energy 
is the mean field dynamo theory [2] . The standard dynamo equation for the large scale magnetic 
field B(i,x) reads <9B/<9i = curl(aB) + /3AB+curl(u x B), where u is the mean velocity field, 
a. is the coefficient of the a-effect and (3 is the turbulent magnetic diffusivity. This equation has 
been widely used for analyzing the generation of the large-scale magnetic field. Traditionally the 
mathematical procedure consists of looking for exponentially growing solutions of the dynamo 
equation with appropriate boundary conditions. While this approach has been quite successful 
in the prediction of large scale magnetic field generation, it fails to predict the subcritical onset 
of a large-scale magnetic field for some turbulent flow. Although the trivial solution B = is 
linearly stable for the subcritical case, the non-normality of the linear operator in the dynamo 
equation for some turbulent flow configurations leads to the transient growth of initial pertur- 
bations It turns out that the non-linear interactions and random fluctuations might amplify 
this transient growth further. Thus, instead of the generation of the large scale magnetic field 
being a consequence of the linear instability of trivial state B = 0, it results from the interaction 
of transient amplifications due to the non-normality with nonlinearities and stochastic pertur- 
bations. The importance of the transient growth of magnetic field for the induction equation 
has been discussed recently in [3 IB]- Comprehensive reviews of subcritical transition in hydro- 
dynamics due to the non-normality of the linearized Navier-Stokes equation, and the resulting 
onset of shear flow turbulence, can be found in |H] ITU]. 

The main purpose of this Letter is to study the non-normal and stochastic amplification of 
the magnetic field in galaxies. Our intention is to discuss the generation of the large-scale mag- 
netic field that cannot be explained by traditional linear eigenvalue analysis. It is known that 
non-normal dynamical systems have an extraordinary sensitivity to stochastic perturbations that 
leads to great amplifications of the average energy of the dynamical system . Although the 
literature discussing the mean field dynamo equation is massive, the effects of non-normality and 
random fluctuations are relatively unexplored. Several attempts have been made to understand 
the role of random fluctuations in magnetic field generation. The motivation was the observation 
of rich variability of large scale magnetic fields in stars and galaxies. Small scale fluctuations 
parameterized by stochastic forcing were the subject of recent research by Farrell and Ioannou 
[7]. They examined the mechanism of stochastic field generation due to the transient growths 
for the induction equation. They did not use the standard closure involving a and (3 param- 
eterization. Hoyon with his colleagues has studied the effect of random alpha-fluctuations on 
the solution of the kinematic mean- field dvnamo|12|. However they did not discuss the non- 
normality of the dynamo equation and the possibility of stochastic transient growth of magnetic 
energy. Both attempts have involved only the linear stochastic theory. Numerical simulations of 
magnetoconvection equations with noise and non-normal transient growth have been performed 
in0 

It is the purpose of this Letter to present a simple stochastic dynamo model for the thin-disk 
axisymmetric aO, dynamo involving three factors: non-normality, non-linearity and stochastic 
perturbations. Recently it has been found JB] that the interactions of these factors leads to 
noise-induced phase transitions in a "toy" model mimicking a laminar-to-turbulent transition. 
In this Letter we discuss three possible mechanisms for the generation of a magnetic field that are 
not based on standard linear eigenvalue analysis of the dynamo equation. The first mechanism 
is a deterministic one that describes an interplay between linear transient growth and nonlinear 
saturation of both turbulent parameters: a and (3. The second and third are stochastic mecha- 
nisms that account for the interaction of the non-normal effect generated by differential rotation 
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with random additive and multiplicative fluctuations. 

Here we study the nonnormality and stochastic perturbation effects on the growth of galactic 
magnetic field by using a Moss's "no-z" model for galaxies Despite its simplicity the "no-z" 
model proves to be very robust and gives reasonable results compared with real observations. We 
consider a thin turbulent disk of conducting fluid of uniform thickness 2h and radius R (R^> h), 
which rotates with angular velocity f2(r) [HIE! We consider the case of afi— dynamo for which 
the differential rotation dominates over the a-effect. Neglecting the radial derivatives one can 
write the stochastic equations for the azimuthal, B v (t) , and radial, B r (t) , components of the 
axisymmetric magnetic field 

dB r a(\B\,Ut)) ^ 2 /?(|B|) 

= h — B * " —w~ Br +if{tl 



dt ~ 9 r 4/i 2 



gB r t ^- u -B ip , (1) 



where a(|B|, is the random non-linear function describing the a— effect, /3(|B|) is the 

turbulent magnetic diffusivity, g = rd£l/dr is the measure of differential rotation (usually 
rdQ/dr < 0). 

Nonlinearity of the functions a(|B|, and /3(|B|) reflects how the growing magnetic field 

B affects the turbulent dynamo coefficients. This nonlinear stage of dynamo theory is a topic 
of great current interest, and, numerical simulations of the non- linear magneto- hydrodynamic 
equations are necessary to understand it. There is an uncertainty about how the dynamo 
coefficients are suppressed by the mean field and current theories seem to disagree about the 
exact form of this suppression [Ej- Here we describe the dynamo saturation by using the 
simplified forms [3] 

a(\B\,Ut) = H + Ut))fa(\B\), /3(|B|) = Po<p p (\B\), (2) 
where (/2 aij g(|B|) is a decaying function such that </? Qj( g(0) = 1. In what follows we use [Sj 

^(|B|) = (l + UB^/B^y 1 , ^(|B|) = (l + 1 + (J ^ /jBy)2 ) , (3) 

where k a and kp are constants of order one, and B eq is the equipartition strength. It should 
be noted that for the afi— dynamo the azimuthal component B v (t) is much larger than the 
radial field B r (t) , therefore, B 2 ~ I? 2 . We did not include the strong dependence of a and f3 on 
the magnetic Reynolds number R m . The back reaction of the magnetic field on the differential 
rotation is also ignored. 

The multiplicative noise £ a (i) describes the effect of rapid random fluctuations of a. We 
assume that they are more important than the random fluctuations of the turbulent magnetic 
diffusivity f3 \l'2\ . The additive noise represents the stochastic forcing of unresolved scales 
[7]. Both noises are independent Gaussian random processes with zero means < £ a (t) >= 0, 
< (,f(t) >= and correlations: 

< UVUs) >= 1D a 5{t - s), < C/(t)e/(«) >= 2D f 5(t - s). (4) 

The intensity of the noises is measured by the parameters D a and D f . One can show ^S] that 
the additive noise in the second equation in (^Q) is less important. 
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The governing equations © can be nondimensionalized by using an equipartition field 
strength B eq , a length h, and a time Oq , where fio is the typical value of angular velocity. 
By using the dimensionless parameters 

n I i x Ra 71-2 d a h Q h 2 . . 

g -> -il \g\, = —, £ = ——, R a = ——, R (jJ = —— : (5) 

/)„ l'i„ p P 

we can write the stochastic dynamo equations in the form of SDE's 

dB r = -(8(p a (B v )B v + e(pp(B v )B r )dt - \f2a[y a (B 9 )B v dWi + V2^dW 2 , 

dB v = -{\g\B r + eip p {B v )B v )dt , (6) 

where W\ and W2 are independent standard Wiener processes. The dynamical system © is 
subjected to the multiplicative and additive noises with the corresponding intensities: 

D a D f 

It is well-known that the presence of noise can dramatically change the properties of a dynamical 
system ^S]. Since the differential rotation dominates over the a-effect (R a <C the system 

© involves two small parameters S = Ra/R^ and e = 1/Ru whose typical values are 0.01 — 0.1 
= 10 — 100, R a = 0.1 — 1). These parameters play very important roles in what follows. 
For small values 5 and e , the linear operator in © is a highly non-normal one ( \g\ ~ 1). This 
can lead to a large transient growth of the azimuthal component B v (t) in a subcritical case. We 
then expect a high sensitivity to stochastic perturbations. Similar deterministic low-dimensional 
models have been proposed to explain the subcritical transition in the Navier-Stokes equations 
(see, for example, |15U16| ). The main difference is that the nonlinear terms in © are not energy 
conserving. 

The probability density function p(t, B r , B v ) obeys the Fokker-Planck equation associated 
with © OH 

dv d d 

9t = ~QB r [{fy(B v )B v +eip{Bv)B r )p] - — \(\g\B r + s^B^B^p] + 

'OBf 

Using this equation in the linear case one can find a closed system of ordinary differential 
equations for the moments < B 2 >, < B r B^ >, and < B 2 > 



(a lV 2 (B v )B 2 +a 2 )-^. (8) 



d_ 

It 




<B 2 > \ I -2e 
< B r B v > \ = \ -\g\ -2e -5 j [ < B r B v > | + | | . (9) 
<B 2 > J { 
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The linear system of equations Q allows us to determine the initial evolution of the average 
magnetic energy E(t) =< B 2 > + < B r B v > + < B 2 > . Similar equations emerge in a variety 
of physical situations, such as models of stochastic parametric instability that explain why the 
linear oscillator subjected to multiplicative noise can be unstable |20| . 

Now we are in a position to discuss three possible scenarios for the subcritical generation of 
galactic magnetic field. 
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Figure 1: Linear case: the azimuthal component B v as a function of time (-B^(O) = 0). for 
\ g \ = 1,6= 1CT 4 and e = 2 • 10~ 2 and different initial values of B r : -0.017, -0.021, -0.03. 

Deterministic subcritical generation. Let us examine the deterministic transient growth 
of the magnetic field in the subcritical case. To illustrate the non-normality effect consider first 
the linear case without noise terms. The dynamical system ((BJ takes the form 




-e 
1.9 1 




(10) 



Since 6 « 1, e « 1 and \g\ ~ 1, this system involves a highly non- normal matrix. Even in 
the subcritical case (0 < 6 < £ 2 /\g\ see below,) when all eigenvalues are negative, B v exhibits 
a large degree of transient growth before the exponential decay. Assuming that B r {t) = e 7 *and 
Btp(0) = b e 7 *we find two eigenvalues 71.2 = — e ± \/<5|g| (the corresponding eigenvectors are 
almost parallel). The supercritical excitation condition 71 > can be written as v^M > £ or 
\/R a Ruj\g\ > vr 2 /4 4 . Consider the subcritical case when < 6 < £ 2 /\g\- The solution of the 
system (fTU|) with the initial conditions B r (0) = — 2c^5/\g\, B^(0) = is 

B r {t) = -cJ ^-(e^ + e^), B v (t) = c(e^ - e^). (11) 



Thus Bp(t) exhibits large transient growth over a timescale of order \je before decaying ex- 
ponentially. In Fig. 1 we plot the azimuthal component B v as a function of time for \g\ = 1, 
6 = 10" 4 and e = 2 • 10~ 2 and different initial values of B r (B^O) = 0). 

Of course without nonlinear terms any initial perturbation decays. However if we take into 
account the back reaction suppressing the effective dissipation (^(|B|) is a decaying function), 
one can expect an entirely different global behaviour. In the deterministic case there can be three 
stationary solutions to © • In Fig. 2 we illustrate the role of transient growth and nonlinearity in 
the transition to a non-trivial state using © with k a = 0.5 and kn = 3. We plot the azimuthal 
component B^ as a function of time with the initial condition -B v (0) = 0. We use the same 
values of parameters \g\, 6 and e and three initial values of B r (0) as in Fig. 1. One can see from 
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Figure 2: Nonlinear case: the azimuthal component B^ as a function of time (B v (0) = 0). for 
\ g \ = 1,6= 10" 4 and e = 2 • 1(T 2 and different initial values of B r : -0.017, -0.021, -0.03. 

Fig. 2 that the trivial solution B v = B r = is nonlinearly unstable to small but finite initial 
perturbations of B r , such as, B r (0) = —0.03. For fixed values of the parameters in nonlinear 
system (J5J, there exists a threshold amplitude for the initial perturbation, above which B v (t) 
grows and below which it eventually decays. 

Stochastic subcritical generation due to additive noise. This scenario has been 
already discussed in the literature jjj (see also for hydrodynamics) . The physical idea is that 
the average magnetic energy is maintained by additive Gaussian random forcing representing 
unresolved scales. It is clear that the non-zero additive noise [a 2 7^ 0) ensures the stationary 
solution to ((§])■ If we assume for simplicity o\ = and 5 = then the dominant stationary 
moment is 

< K >- £?■ m 

We can see that due to the non-normality of the system (|10|) the average stationary magnetic 
energy E s t ~< B^ > s t exhibits a high degree of sensitivity with respect to the small parameter 
e:E st ~ £ ~ 3 [TU Ell- 
Stochastic subcritical generation due to multiplicative noise. Here we discuss the 
divergence of the average magnetic energy E(t) =< B^. > + < BrB^ > + < B^ > with time t 
due to the random fluctuations of the a— parameter. Although the first moments tend to zero in 
the subcritical case, the average energy E(t) grows as e xt when the intensity of noise a\ exceeds 
a critical value. The growth rate A is the positive real root of the characteristic equation for the 
system © 

(A + 2ef - 4%|(A + 2e) - 2<n\g\ = 0. (13) 

For 6 = 0, the growth rate is Ao = — 2e + (2<ri|g|) 1 / 3 as long as it is positive, and the excitation 
condition can be written as a± > a cr = 4e 3 /\g\. It means that the generation of average magnetic 
energy occurs for ao = ! It is interesting to compare this criterion with the classical supercritical 
excitation condition: 5\g\ > e 2 4 . To assess the significance of this parametric instability it is 
useful to estimate the magnitude of the critical noise intensity a cr . First let us estimate the 
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parameter e = ir 2 f3 / (AQoh 2 ) . The turbulent magnetic diffusivity is given by /? ~ lv/3, where v 
is the typical velocity of turbulent eddy v — 10 km s _1 , and I is the turbulent scale, I ~ 100 
pc. For spiral galaxies, the typical values of the thickness, h, and the angular velocity, Qq, are 
h ~ 800 pc and S7o — 10~ 15 s -1 ; \g\ ~ 1 4 j. It gives an estimate for e ~ 3.2 x 10~ 2 , that 
is, a cr ~ 1.3 x 10~ 4 . In general X(5) = A + (4/3|c/|(2cti \g\)~ 1/3 )5+ o(S). This analysis predicts 
an amplification of the average magnetic energy in a system (JSJ) where no such amplification is 
observed in the absence of noise. The value of the critical noise intensity parameter a cr , above 
which the instability occurs, is proportional to e 3 , that is, very small indeed. To some extent, the 
amplification process exhibits features similar to those observed in the linear oscillator submitted 
to parametric noise j2H] • To avoid the divergence of the average magnetic energy, it is necessary 
to go beyond the kinematic regime and consider the effect of nonlinear saturations. 

In summary, we have discussed galactic magnetic field generation that cannot be explained by 
traditional linear eigenvalue analysis of dynamo equation. We have presented a simple stochastic 
model for the a£l dynamo involving three factors: (a) non-normality due to differential rota- 
tion, (b) nonlinearity of the turbulent dynamo a effect and diffusivity /?, and (c) additive and 
multiplicative noises. We have shown that even for the subcritical case, there are three possible 
scenarios for the generation of large scale magnetic field. The first mechanism is a determin- 
istic one that describes an interplay between transient growth and nonlinear saturation of the 
turbulent a— effect and diffusivity. We have shown that the trivial state B = can be non- 
linear ly unstable with respect to small but finite initial perturbations. The second and third 
are stochastic mechanisms that account for the interaction of non-normal effect generated by 
differential rotation with random additive and multiplicative fluctuations. We have shown that 
multiplicative noise associated with the a— effect leads to exponential growth of the average 
magnetic energy even in the subcritical case. 
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